Memory for Light as a Quantum Process 
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We report complete characterization of an optical memory based on electromagnetically induced 
transparency. We recover the superoperator associated with the memory, under two different work- 
ing conditions, by means of a quantum process tomography technique that involves storage of 
coherent states and their characterization upon retrieval. In this way, we can predict the quantum 
state retrieved from the memory for any input, for example, the squeezed vacuum or the Fock state. 
We employ the acquired superoperator to verify the nonclassicality benchmark for the storage of a 
Gaussian distributed set of coherent states. 

PACS numbers: 42.50.Ex, 03.67.-a, 32.80.Qk, 42.50.Dv 



Introduction Quantum memory for light is an essen- 
tial technology for long distance quantum communica- 
tion [l[ and for any future optical quantum informa- 
tion processor. Recently, several experiments have shown 
the possibility to store and retrieve nonclassical states of 
light such as the single photon 0, 0] , entangled [1] and 
squeezed vacuum Q states using coherent interactions 
with an atomic ensemble. 

In order to evaluate the applicability of a quantum 
memory apparatus for practical quantum communication 
and computation, it is insufficient to know its perfor- 
mance for specific, however complex, optical states, be- 
cause in different protocols, different optical states are 
used for encoding quantum information 0, 0] ■ Practical 
applications of memory require answering a more general 
question: how will an arbitrary quantum state of light be 
preserved after storage in a memory apparatus? 

Here we answer this question by performing complete 
characterization of the quantum process associated with 
optical memory based on electromagnetic induced trans- 
parency (EIT) 0. Memory characterization is achieved 
by storing coherent states (i. e. highly attenuated laser 
pulses) of different amplitudes and subsequently measur- 
ing the quantum states of the retrieved pulses. Based 
on the acquired information, the retrieved state for any 
arbitrary input can be predicted and additionally, any 
theoretical benchmark on quantum memory performance 
can be readily verified. 

Coherent state quantum process tomography We can 
define complete characterization of an optical quantum 
memory as the ability to predict the retrieved quantum 
state £ (p) when the stored input state p is known. This 
is a particular case of the quantum "black box" problem, 
which is approached through a procedure called quan- 
tum process tomography (QPT) [§]. QPT is based on 
the fact that every quantum process (in our case, opti- 
cal memory) is a linear map on the linear space L(H) of 
density matrices over the Hilbert space H on which the 
process is defined. The associated process can thus be 
characterized by constructing a spanning set of "probe" 
states in L(H) and subjecting each of them to the action 



of the quantum "black box" . If we measure the process 
output £{p{) for each member pi of this spanning set, 
we can calculate the process output for any other state 
p = J2i a iPi according to 



(1) 



The challenge associated with this approach is the con- 
struction of the appropriate spanning set, given the infi- 
nite dimension of the optical Hilbert space and the lack 
of techniques for universal optical state preparation. For 
this reason, characterizing memory for light, that is not 
limited to the qubit subspace, is much more difficult than 
memory for superconducting qubits, which has been re- 
ported recently [13] ■ Our group has recently developed a 
process characterization technique that overcomes these 
challenges . Any density matrix p of a quantum opti- 
cal state can be written as a linear combination of density 
matrices of coherent states \a) according to the optical 
equivalence theorem 



p = 2 J Pp(a)\a){a\d 2 a, 



(2) 



where Pp(a) is the state's Glauber-Sudarshan P-function 
and the integration is performed over the entire com- 
plex plane. Although the P-function is generally highly 
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FIG. 1: (color online). Schematic of the experimental setup 
used to characterize the process associated with the quantum 
memory. PBS, polarizing beam splitter. 
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singular, any quantum state can be arbitrarily well ap- 
proximated by a state with an infinitely smooth, rapidly 
decreasing P-function [T2|. Therefore, by measuring how 
the process affects coherent states, one can predict its 
effect on any other state. The advantage of such ap- 
proach (which we call coherent-state quantum process 
tomography or csQPT) is that it permits complete pro- 
cess reconstruction using a set of "probe" states that are 
readily available from a laser. 

Experimental setup We performed csQPT on opti- 
cal memory Q realized in a warm rubidium vapor 
by means of electromagnetically-induced transparency 
(Fig. [T]). The atoms are 87 Rb and the vapor tempera- 
ture is kept constant at 65° C. 

The signal field is resonant with the 5 S1/2, F = 1) <-> 
\ 5 Pi/2,F = 1) transition at 795 nm and is produced by 
a continuous-wave Ti:Sapphire laser. An external cavity 
diode laser, phase locked at 6834.68 MHz to the signal 
laser 13| serves as the EIT control field source, and is 



resonant with the \ 5 Si/ 2 ,F = 2) «-> \ 5 Pi/ 2 ,F = 1) tran- 
sition. The fields are red detuned from resonance by 630 
MHz in order to improve the storage efficiency. The con- 
trol field power is 5 mW and the beam spatial profile is 
mode matched with the signal beam to a waist of 0.6 mm 
inside the rubidium cell. Signal and control fields are or- 
thogonally polarized; they are mixed and separated using 
polarizing beam splitters. 

The two photon detuning A2 between the signal and 
control fields is modified by varying the frequency of the 
control field laser through the phase lock circuit, while 
an acousto-optical modulator (AOM) is used to switch 
on and off the control field intensity. We analyzed two 
different operative conditions characterized by A2 = 
and 0.54 MHz. 

The input pulse is obtained by chopping the 
continuous-wave signal beam via an AOM to produce 
1 /is pulses [Fig. [2jc)] with a 100 kHz repetition rate. 
A second AOM is used to compensate for the frequency 
shift generated by the first. Transfer of the light state 
into the atomic ground state superposition (atomic spin 
wave) is accomplished by switching the control field off 
for the storage duration of r = 1 ^s when the input pulse 
is inside the rubidium cell. 

We performed full state reconstruction of both the in- 
put and retrieved fields by time domain homodyne to- 
mography A part of the TkSapphire laser beam 
serves as a local oscillator for homodyne detection; while 
its phase is scanned via a piezoelectric transducer, the 
homodyne current is recorded with an oscilloscope. For 
every state, 50000 samples of phase and quadrature are 
measured and processed by the maximum likelihood al- 
gorithm TH, 3|, estimating the state density matrix in 
the Fock basis. 

Tomography of quantum memory In order to deter- 
mine the coherent state mapping necessary for recon- 
structing the process, we measured 10 different coherent 
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FIG. 2: (color online). Wigner functions of input coherent 
states with two different amplitudes and the corresponding 
retrieved states £(|ai)(ai|) and £ (\az)((X2\). The input state 
amplitudes are a± = 2.3, Q2 = 10.3 (a) and ai — 4.5, 
02 = 10.9 (b). Two-photon detunings between the control 
and signal fields are 0.54 MHz (a) and MHz (b). Input 
pulse (black dashed line) and retrieved light (red solid line) 
(c). Variance of the phase quadrature as a function of the 
retrieved state amplitude (d). 



states \oti) with mean photon numbers ranging from 
to 285 along with their corresponding retrieved states 
£{\ai){oti\) [Fig. [2ja) and (b)]. Subsequently, we ap- 
plied polynomial interpolation to determine the value of 
£{\a){a\) for any value of a in the range to 16.9. Per- 
forming tomographic reconstruction for these highly dis- 
placed states requires good phase stability between the 
signal and local oscillator. Phase fluctuations produce an 
artefact in the reconstruction in the form of amplitude 
dependent increase in the phase quadrature variance. In 
our measurements, the reconstructed input states |a,-) re- 
semble theoretical coherent states with a fidelity higher 
than 0.999 for mean photon values up to 150 [Fig. [21a) 
and (b)]. 

By inspecting the Wigner functions of the input and 
retrieved states, one can clearly notice the detrimental 
effects of the memory. First, there is attenuation of the 
amplitude by a factor of 0.41 ± 0.01 for the signal field in 
two-photon resonance with the control, which increases 
to a factor of 0.33 ± 0.02 when a two-photon detuning 
of A2 = 0.54 MHz is introduced. This corresponds to a 
mean photon number attenuation by factors of 0.17±0.02 
and 0.09 ± 0.01, respectively. Note that in the case of 
nonzero two-photon detuning, the attenuation is greater 
than the factor of 0.14 obtained in classical intensity mea- 
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FIG. 3: (color online). The diagonal elements of the process 
tensor Xfcfc"\ measured by csQPT in the Fock basis for A2 
(a) and 0.54 MHz (b). 



surcment [Fig. [2c)]. This is because the temporal mode 
of the retrieved state is slightly chirped, and could not 
be perfectly matched to the mode of the local oscillator. 

Second, retrieved coherent states experience an in- 
crease in the phase quadrature variance that depends 
quadratically on the state amplitude. This effect pro- 
duces an ellipticity in the retrieved state Wigner function 
(Fig. [^a) and (b)] and can be attributed to the noise in 
the phase lock between the signal and control lasers [l3| • 
Fluctuations A<fi of the relative phases between the two 
interacting fields randomize the phase of the retrieved 
signal field with respect to the local oscillator. Assum- 
ing a Gaussian distribution for A(f> with zero mean and 
variance the variance of the phase quadrature can be 
expressed as: 



a « = 2 + T ' 1 



e 2 < 



(3) 



where qg is the mean amplitude. We fit our experimental 
data with Eqj3]and estimate an II standard deviation 
for A(j> [Fig. HTd)1, in agreement with independent esti- 
mates fl3| . 

The third detrimental effect preventing the atomic en- 
semble from behaving as a perfect memory is the pop- 
ulation exchange between atomic ground states (ItI Il8|. 
Besides limiting the memory lifetime, this exchange gen- 
erates spontaneously emitted photons in the signal field 
mode adding an extra noise that thermalizes the stored 
light by increasing the quadrature variance indepen- 
dently of the input amplitude and phase. We measured 
the extra noise from the quadrature variance of retrieved 
vacuum states £(|0)(0|) and found it to equal 0.185 dB 
when both fields were tuned exactly at the two photon 
resonance, which corresponds to the mean photon num- 



Tr 



n£(|0)(0| 



ber in the retrieved mode equal to n 

0.022. This noise is reduced to 0.05 dB (corresponding 
to n = 0.005 ) in the presence of two photon detuning. 
For this reason, it is beneficial to implement storage of 
squeezed light in the presence of two-photon detuning, in 
spite of higher losses. 

In the presence of the two-photon detuning, the evo- 
lution of the atomic ground state superposition brings 



about a phase shift of the retrieved state with respect to 
the input by 27rA 2 r = 200° as is visible in Fig. [Ha). 

Based on the information collected from the storage of 
coherent states, we reconstruct the memory process in 

3' 



the %-matrix representation, defined by 



£{p) = 



E 



n.in 

X k i 



Al, n pA 



rn.ki 



(4) 



where xl '1 ^ s the rank 4 tensor comprising full informa- 
tion about the process and Aij is a set of operators that 
form a basis in the space of operators on EL Since H is 
the Hilbert space associated with an electromagnetic os- 
cillator, it is convenient to choose A^j = where \i) 
and \ j) are the photon number states. The details of cal- 
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culating the process tensor are described elsewhere 
Fig. [3] displays the diagonal subset xT ™ °f the P rocess 
tensor elements. 

Performance tests In order to verify the accuracy of 
our process reconstruction, we have used it to calculate 
the effect of storage on squeezed vacuum with A2 = 0.54 
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FIG. 4: (color online). Comparison of the experimentally 
measured squeezed vacuum states retrieved from the quantum 
memory and those predicted with csQPT. For each case, the 
Wigner function and the quadrature variance as a function 
of the local oscillator phase are shown, (a), Experimental 
measurement [f| with A2 = 0.54 MHz. (b), Prediction with 
A2 = 0.54 MHz. (c), Experimental measurement with A2 = 
0. (d), Prediction with A 2 = 0. 
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MHz, as studied in a recent experiment of our group [f|, 
and with A 2 = MHz. We applied the superoperator 
tensor measured with csQPT to the squeezed vacuum 
produced by a subthreshold optical parametric ampli- 
fier with a noise reduction in the squeezed quadrature 
of —1.86 dB and noise amplification in the orthogonal 
quadrature of 5.38 dB (i.e. the same state as used as the 
memory input in Ref. Q). In this way, we obtained a 
prediction for the state retrieved from the memory, which 
we then compared with the results of direct experiments. 
This comparison yields quantum mechanical fidelities of 
0.9959 ± 0.0002 and 0.9929 ± 0.0002 for the two-photon 
dctunings of A 2 = 0.54 MHz and A 2 = respectively 
(Fig.©. 

As discussed above, zero detuning warrants lower 
losses (thus higher amplitude of the noise variance) and 
no phase rotation, but higher excess noise (thus no 
squeezing in the retrieved state). Nevertheless the two 
photon resonant configuration offers a better fidelity if 
the single photon state is stored 0, S| , as evidenced by 
comparing the superoperator element x\'i 01 Fig. [3] (a) 
and (b). 

In addition to the ability to predict the output of the 
memory for any input state, our procedure can be used 
to estimate the performance of the memory according to 
any available benchmark. As an example, we analyze the 
performance of our memory with respect to the classical 
limit on average fidelity associated with the storage of co- 
herent states with amplitudes distributed in phase space 
according to a Gaussian function of width 1/A [2l[. This 
limit as a function of A is given by: 

H-oo 1 + A 

F(X) = 2A / cxp (-Xa 2 )(a\£ (\a) (a\)\a) ada < 

Jo 2 + A 

(5) 

From csQPT data, we evaluate the average fidelity asso- 
ciated with our memory for both values of A 2 (Fig. O. 
Both configurations show nonclassical behavior. The 
higher value of average fidelity correspond to A 2 = 
and is explained by a higher storage efficiency. 

Conclusion In summary, we have demonstrated com- 
plete characterization of an EIT-bascd quantum memory 
by csQPT. This procedure allows one to predict the effect 
of the memory on an arbitrary quantum-optical state, 
and thus provides the "specification sheet" of quantum- 
memory devices for future applications in quantum in- 
formation technology. Furthermore, our results offer in- 
sights into the detrimental effects that affect the storage 
performance and provide important feedback for the de- 
vice optimization. We anticipate this procedure to be- 
come standard in evaluating the suitability of a memory 
apparatus for practical quantum telecommunication net- 
works. 
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FIG. 5: (color online). Average fidelity of the quantum mem- 
ory for a Gaussian distributed set of coherent states. Blue 
empty (red filled) dots show the average fidelity calculated 
from the csQPT experimental data for A2 = (0.54 MHz). 
The experimental uncertainty is 0.0002. The solid line shows 
the classical limit [2l| . 
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